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In light of the discovery of a Higgs-like particle at the LHC, we revisit the status of the precision 
electroweak data, focusing on two discrepant observables: 1) the long-standing 2.4a deviation in 
the forward-backward asymmetry of the bottom quark A b FB , and 2) the 2.3a deviation in Rb, the 
ratio of the Z — > bb partial width to the inclusive hadronic width, which is now in tension after a 
recent calculation including new two-loop electroweak corrections. We consider possible resolutions 
of these discrepancies. Taking the data at face value, the most compelling scenario is that new 
physics directly affects A b FB and Rb, bringing the prediction into accord with the measured values. 
We propose a modified 'Beautiful Mirrors' scenario which contains new vector-like quarks that mix 
with the b quark, modifying the Zbb vertex and thus correcting A b FB and Rb- We show that this 
scenario can lead to modifications to the production rates of the Higgs boson in certain channels, and 
in particular a sizable enhancement in the diphoton channel. We also describe additional collider 
tests of this scenario. 



The ATLAS and CMS experiments have recently dis- 
covered a new boson with properties closely resembling 
those of the Standard Model (SM) Higgs boson [TJ [S]. 
The focus of the experimental collaborations now turns 
to determining the properties of this state, such as its 
quantum numbers and its couplings to other SM parti- 
cles. This is an extremely important task since any de- 
viation from the SM Higgs boson predictions will unam- 
biguously point to new physics (NP) at the TeV scale. In 
this respect it is intriguing that both experiments observe 
a slight enhancement in the h — ¥ 77 channel, though with 
the current dataset this enhancement is not statistically 
significant [21 0]. 

If, as the current data suggest, this new state is in- 
deed the Higgs boson, its mass rrih ~ 125 GeV will be 
in accord with the expectation indirectly suggested by 
the precision electroweak data, which favors a light SM- 
like Higgs boson. Of the many precision measurements 
used to test the electroweak sector, most of them are 
in good agreement with the SM predictions, see for ex- 
ample [5 , 6 ,[7 . However, there are a couple of notable 
discrepancies. The forward-backward asymmetry of the 
bottom quark A b FB measured at the Z-pole at LEP1 con- 
stitutes a long-standing exception. The measured value 
and SM prediction are [5], [7], 

(^fb)cxp = 0.0992 ±0.0016, (1) 
(^fb)sm = 0.1032+°;ooog, (2) 

and thus exhibit a 2.4er discrepancy. Furthermore, a re- 
cent calculation of Rb, the ratio of the Z — > bb partial 
width to the inclusive hadronic width, which includes 
new two-loop electroweak corrections, now puts the pre- 
diction in tension with the measured value [5] . The mea- 
surement and prediction read [5], [7], 

(R b ) cxp = 0.21629 ± 0.00066, (3) 
(-Rh)sM = 0.21474 ± 0.00003, (4) 



displaying a 2.3<r discrepancy. 

It is a matter of debate whether these deviations call 
for NP. Among a large ensemble of measurements, such 
as that used in the electroweak fit, one may expect an 
occasional discrepancy of this size, which could be the 
result of a statistical fluctuation or unknown systematic 
error. However, a possible counterpoint with regards to 
A b FB , as emphasized in Refs. [91 [TU], is that A b FB fa- 
vors a heavy Higgs. Excluding the measurement from 
the global fit would favor a Higgs boson much lighter 
than the LEP bound, perhaps suggesting NP that im- 
proves the electroweak fit with a 125 GeV Higgs. This 
approach was also advocated in Ref. [TT]. This dis- 
crepancy and its interplay with the indirect determi- 
nation of the Higgs mass has stimulated a number of 
theoretical works that attempt to resolve this mystery, 

see e.g. mmmmmmmmmmm, 
mm- 

With a direct and precise measurement of the Higgs 
mass from the LHC experiments, as well as the relatively 
recent tension in Rb, it is useful to revisit these argu- 
ments. At a very simplistic level, we observe that there 
are now two measurements related to bottom quarks 
which are discrepant, suggesting either a systematic ef- 
fect associated with 6-quarks or possibly NP in the b- 
quark sector. Thus, we will consider two possible reso- 
lutions to the A b Fg and Rb discrepancies: First, we will 
entertain the possibility that NP alters the Zbb vertex, 
thus changing the predictions for A b FB and Rb- As we 
will demonstrate, this scenario leads to a dramatic im- 
provement in the global fit of the data compared to the 
SM. Second, we will hypothesize that the A b FB and Rb 
deviations are a result of an unknown systematic effect, 
and as such remove these observables from the EW fit. In 
this case, we find that the remaining observables are well 
described by the SM, and NP which contribute to oblique 
parameters S and T |24j can only marginally improve the 



2 



fit. 

Thus, if one takes the A b FB and Rb measurements se- 
riously, the most compelling interpretation is that NP 
directly affects A b FB and R b so as to bring the measured 
and predicted values into agreement. As an example we 
will focus on the "Beautiful Mirrors" scenario of Ref. [13] • 
In this class of models, new vector-like quarks mix with 
the bottom quark of the SM, thus modifying the Zb R b R 
coupling and in turn A b FB . We then expect that the Higgs 
production rates as well as decay modes will deviate from 
their SM values, since the hbb coupling is altered and new 
exotic vector-like quarks contribute to the loop induced 
processes h — > gg and h — > 77. We will revisit minimal 
models of this type, as well as some extensions to exhibit 
the range of possibilities in Higgs phenomenology. 

INTERPRETING THE PRECISION 
ELECTROWEAK DATA 

In this section we will examine interpretations of the 
precision electroweak data using global fits, particularly 
focusing on possible solutions to the A h FB and R b dis- 
crepancies. Specifically, there are two distinct scenarios 
we wish to test: 

• New physics affects the Zbb vertex, thus modifying 
the A b FB and Rb- 

• The discrepancies in A b FB and Rb are attributed to 
an unknown systematic error, and the observables 
are thus removed from the global fit. 

Before presenting our results for these two scenarios, let 
us briefly describe the ingredients going into our global fit 
of the precision electroweak data. Our fit follows closely 
that of the Gfitter group [7J. We use the high energy 
data, including Z-pole data from LEP and SLD [5], as 
well as measurements of the top|2,5j , W and Higgs 
mass [I] , [2] and the hadronic contribution to the running 
fine structure constant |27j . Explicitly, we consider the 
following inputs in the fit: 

mz, Tz, cr£ ad , Re, R c , Rb, 
A F B , A(, A c , Ab, A FB , A b FB , sin 2 9 e g, 

m w , IV, m t , Aa^) m h- (5) 

The experimental values for these observables are taken 
from [7J. The leptonic asymmetry parameter An is an 
average of the LEP and SLD measurements, while the 
heavy flavor observables A c , Ab, A c FBl A b FB , R c , Rb 
are averaged from LEP and SLC[S]. We also account for 
correlations in Z lineshape as well as heavy flavor observ- 
ables [5]. The W boson mass is taken from [26], which 
includes the recent precise measurements from CDF and 
DO, as well as older measurements from LEP and the 
Tevatron experiments. For the Higgs mass, we combine 



the ATLAS (126.0 ± 0.4 (stat) ± 0.4 (syst) GeV pQ) and 
CMS (125.3 ± 0.4 (stat) ± 0.5 (syst) GeV 0) results, as- 
suming no correlation of the systematic uncertainties, to 
obtain m h = 125.7 ± 0.4 GeV. 

For the SM theory predictions we allow the following 
parameters to vary in the fit: 

m h ,m z ,m t , Aaf^,a s . (6) 

We work with the running electromagnetic coupling 
which depends on the Aa[ 5 a ' d in For the W boson 
mass prediction we utilize the numerical parameteriza- 
tion in Eqs. (6,7,9) of Ref. |28j . which includes the full 
two loop prediction plus certain important higher order 
corrections. We employ the numerical parameterizations 
of the Z boson partial widths in terms of the SM pa- 
rameters ([6| from Ref. [29j (updated from earlier work 
in Refs. [SO]). From here we construct T z , c° ad , Re, R c - 
For Rb , we instead use the recent result of Ref. [8 , which 
contains new two-loop electroweak corrections. For the 
effective weak mixing angle, we use the parameterizations 
in Refs. 131] . which account for the latest two- loop 
electroweak corrections. From these parameterizations, 
we calculate the asymmetry parameters Af and and the 
forward-backward asymmetries A FB . 

Additionally, we will explore NP scenarios with oblique 
corrections parameterized with the Peskin-Takeuchi pa- 
rameters S and T, as well as tree-level shifts to the Zbb 
vertex. In the presence of oblique corrections the theory 
predictions for the observables are modified to have the 
form O = Osm + csS + ctT, where the coefficients q de- 
pend on the input parameters |5|. These shifts have been 
computed in e.g. Ref. [32] . and we employ the expres- 
sions from this reference. We note that for the models 
considered in this paper the contribution to the U pa- 
rameter is small, so we will not consider it in our fits 
here. 

We parameterize the shifts in the Zbb coupling as 

C = ^-Z^[{g Lb + Sg Lb )P L + (g Rb + 6g Rb )P R }b 1 (7) 
cw 

where cw = cos 9w, with 9w the weak mixing angle, 
gLb and g R b are the SM left and right couplings given by 
g Lb = 1/2 + s 2 w /3 ~ -0.42, and g Rb = s 2 w /S ~ -0.077. 
The shifts 5gLb,m represent a non-universal contribu- 
tion due to NP. These shifts enter into the predictions 
for observables in ^ depending on the Z — ► bb partial 
width as well as the asymmetry parameters A b FB and Ab- 
Similar fits to Zbb couplings have been performed in the 
past [33], [33], [HI, [TTJ. 

We will now use the results of global fits to the data to 
examine the two scenarios described at the beginning of 
this section, namely 1) NP in Zbb couplings affects A b FB 
and Rb, and 2) A b FB and Rb are a result of an unknown 
systematic error. 
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Fit 


rrih (GeV) 


PullL4i n l 


PulURJ 




no 


v 2 /d o f 


p- value 


SM (All data) 


125.7 ±0.4 


2.6 


-2.4 






20.5/13 


0.08 


SM + Sg Lb>Rb (+) (All data) 


125.7 ±0.4 


0.7 


0.2 


0.001 ±0.001 


0.016 ±0.005 


9.6/11 


0.57 


SM + 8g Lb , Rb (-) (All data) 


125.7 ±0.4 


0.8 


0.2 


0.001 ±0.001 


-0.170 ±0.005 


9.3/11 


0.60 



TABLE I. New physics in Zbb: Results of the global fits to the precision electroweak observables in Eq. |5J) are displayed for 
1) the SM hypothesis and 2) the hypothesis of NP in Zbb couplings as in Eq. |7b. We show the Higgs mass prediction, the pull 
values for A h FB and R b , and the goodness-of-fit in terms of x 2 /d.o.f. and p-value. For the second hypothesis, we also show the 
best-fit regions for the coupling shifts {8gm, Sgm,). Note that there are two best-fit regions regions (labeled by "±" and "— "). 



New physics in Zbb couplings Here we take se- 
riously the measurements of A b FB and Rb and include 
them in the global fit along with the other observables 
in Eq. Before considering NP, we would first like to 
understand how well the SM describes the data. The re- 
sults of our fit to the SM are presented in Table [I] For the 
best-fit point, we obtain x 2 /d.o.f. = 20.5/13 correspond- 
ing to a p-value is 0.08. The goodness-of-fit is marginal 
due to the discrepancies in A b FB (2.6a) and Rb (2.4c). 
The Higgs mass is nailed down by the LHC measure- 
ment, with a value of rrih = 125.7 ± 0.4 GeV. 

We now hypothesize the existence of NP that causes 
non- universal shifts in the Zbb vertex as in Eq. Q. The 
results of the fit are displayed in Table |TJ The quality 
of the fit is much improved compared to the SM fit. 
There are two preferred regions in the fit, one charac- 
terized by a large negative shift in the right coupling 
Sgm = — 0.17 ± 0.05 (labeled "-") and the other charac- 
terized by a moderate positive shift of the right coupling 
Sgm = 0.016 ± 0.005 (labeled "+")■ In both cases, the 
data also suggest a small positive shift in the left coupling 
Sgib = 0.001 ±0.001. The improvement in the fit quality 
stems from the fact that the A b FB and Rb predictions are 
no longer discrepant, lying within la of their measured 
values. In Fig. 0, we show the preferred region in the 
SgLbSgm parameter space for the solution with a moder- 
ate positive Sgm shift, which will be relevant below when 
we investigate explicit models of NP that cause shifts to 
the Zbb vertex. 

The asymmetry parameter At, was measured at SLD 
to be (Ah) cxp = 0.923 ± 0.020, which is in agreement 
with our SM prediction (A b ) SM = 0.9346 ± 0.0003, is 
modified by the shifts SgLb,m- For example, at the 
best fit point (Sgm = 0.016, SgLb = 0.001), we predict 
(A b )th = 0.907 ± 0.009, still within la of the measured 
value. The larger error in the theory prediction is mainly 
due to the uncertainty in Sgm- 

A b FB and Rb due to a systematic error Next, we 
consider the possibility that A b FB and Rb are discrepant 
because of some unknown systematic error, perhaps re- 
lated to the measurement and identification of bottom 
quarks. We thus exclude these observables from the fit. 
The results are shown in Table ITT1 The SM describes well 
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FIG. 1. Best-fit region for the Zbb coupling shifts Sg^t,, 8g Rb . 
We show here the region with a moderate positive shift to 
Sgm- The SM corresponds to the green point. 



the rest of the data, with a p-value of 0.67. The direct 
measurement of the Higgs mass by the LHC experiments 
entirely determines the Higgs mass to be mt = 125.7 
GeV, so given the high quality of the fit, it would seem 
that there is no need for NP. 

To see the tendency of the precision electroweak data 
(excluding A b FB and Rb to favor very light Higgs masses, 
we have also performed the fit while excluding the LHC 
measurement of rrih- One observes that the fit selects 
rrih = 70±30 GeV. One interpretation is that the indirect 
determination displays a slight tension with the direct 
determination, and one could then speculate about NP. 

For example, one way to 'reconcile' a heavy Higgs with 
precision electroweak data is to invoke NP in the form of 
oblique parameters S, T. We have thus done a fit to the 
data allowing S and T to vary. We find that the quality 
of the fit is slightly improved compared to the SM case, 
with a p-value of 0.78. 

Summary There is no known reason to attribute 
the A b FB ,Rb discrepancies to a systematic effect, though 
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Fit 


m h (GeV) 


5 


T 


X 2 /d.o.f 


p- value 


SM (No A b FB ,R b ) 


125.7 ±0.4 






9.3/12 


0.67 


SM (No A b FB ,R b ,m h ) 


70 ±30 






5.6/12 


0.92 


SM + S,T (No A b FB ,R b ) 


125.7 ±0.4 


-0.08 ±0.10 


0.0 ± 0.08 


5.6/9 


0.78 



TABLE II. A b FB and R b due to a systematic error: Results of the global fits to the precision electroweak observables in Eq. fl5), 
excluding those observables listed as "No". We consider 1) the SM hypothesis and 2) the hypothesis of NP in the form of 
oblique parameters S,T. We show the Higgs mass prediction and the goodness-of-fit in terms of x 2 /d.o.f. and p-value. For the 
second hypothesis, we also show the best-fit regions for S, T parameters. 



since both measurements rely on b quark identification, it 
is conceivable that there is an unknown systematic effect. 
In any case, removing these discrepant observables from 
the global fit, the remaining data are in excellent agree- 
ment with the SM (quantified by a sizable p- value in the 
corresponding fit), so it is difficult to make a strong case 
for NP. There is only a modest improvement to the fit if 
one invokes the oblique parameters S and T. 

On the other hand, if one takes the data at face 
value, then there exists a preference for some NP which 
gives a non-universal correction to the Zbb couplings. 
In our view this is the most compelling explanation of 
the A b FB , R b discrepancies. Below we will investigate NP 
models which can cause such non-universal shifts. 



CORRELATING A b FB , R b WITH HIGGS DATA 

The models that we will consider to address the dis- 
crepancies in the precision electroweak data, namely A b FB 
and R b , will also cause modifications to the properties of 
the Higgs boson. There are three effects that occur: 1) 
the coupling of the Higgs boson to bottom quarks, yhbb is 
modified, 2) new colored particles contribute to the pro- 
cess h — > gg, and 3) new charged particles contribute to 
the process h — > 77. These three effects can be described 
in a model independent fashion by the quantities, r 7 , r g , 
rb, defined as 



n 



r(ft -» 77) 

T(h -> 77)sm ' 
r(ft -> gg) 

-> 99)sm' 
T(h -» bb) 

T{h -> bb) 



SM 



(8) 
(9) 
(10) 



The experimental collaborations have presented measure- 
ments of the best-fit signal strength /j, of the new state 
in the channel i, defined as 



M = 



a(pp 







o-(pp — > h — > i) 



SM 



(11) 



This dataset can then be used to constrain or test for 
the presence of new physics, as has been done first in 



Ref. [35], after the initial December 2012 announcement 
of an excess at 125 GeV. For recent fits done after the 
discovery of the 125 GeV boson, see Refs. [3B"] . 

In terms of these quantities, the Higgs signal strengths 
are predicted to be 



M77 

fiVV 
^bb 



r g r ~t 


1 - 


hBr b (r 6 


-1) + 


Br 9 (r s - 1) 






r 9 




1 - 


HBr b (r b 


-1) + 


Br s (r s - 1) 






n 




1 - 


HBr b (r fc 


-1) + 


Br ff (r 9 - 1) 



(12) 
(13) 
(14) 



where Br^ denotes the SM Higgs branching ratio in the 
channel i and the signal strengths are defined as 



<j(pp — > h — > i) 
cr(pp -> h -> i)sm 



(15) 



We have done a simple fit to the combined signal strength 
data using the 77, ZZ, WW, and bb channels from AT- 
LAS [J, CMS [T], and the 77, WW, and bb channels from 
the Tevatron experiments [37] . We do not include the tt 
channel in the results presented here, but have checked 
that this does not qualitatively affect our conclusions. 

The best fit point yields (r 7 ,r g ,r b ) = (2.1,1.4,2.1), 
with a x 2 /d-o.f = 6.3/8 versus x 2 / d - oi = H.l/H for 
the SM. However, there is a strong degeneracy in the 
parameters r g and r b , in that similar fit qualities can 
be achieved for different combinations of these parame- 
ters, as is illustrated in Fig. |2]). For example, the points 
(r 7 ,r g ,r b ) = (2.1,0.7,0.7), (2.1,0.8,1), (2.1,1,1.5) give 
similar fit qualities as the best-fit point. One can see this 
very easily in the regime where r bi9 are large, in which 
case the signal strength predictions in Eqs. ( 12|13|14 | be- 
come dependent only on the ratio r b jr g rather than each 
individual parameter r g and r b . This degeneracy can 
eventually be broken by separating the diphoton data 
into categories sensitive to the production mechanism of 
the Higgs boson, in particular the VBF production mode. 
We have checked, however, that the degeneracy is not 
broken with the current dataset sensitive to VBF pro- 
duction due to the limited precision in these channels. 
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FIG. 2. Best-fit regions of the Higgs signal strength data 
(1, 2, 3a) in the r g — rb plane. Here we have marginalized over 
r 7 , and contours of constant r 1 are represented by the red 
lines. There is a quasi- flat direction in the x 2 function along 
Tb ~ 2 r g for large value of rt.g. 



MODELS FOR A b FB AND R b AND HIGGS 
PHYSICS 

We now investigate NP models that address the A b FB 
and R b discrepancies in the precision electroweak data. 
We will consider the "Beautiful Mirrors" scenario [14] as 
our canonical example of NP that directly alters A FB and 
Rb and simultaneously leads to interesting deviations in 
Higgs phenomenology. The basic mechanism for correct- 
ing A b FB and Rb is mixing between the bottom quark of 
the SM and new vector-like 'mirror' quarks. This causes 
a modification of the coupling of the Z boson to the phys- 
ical bottom quark, thus affecting A b FB and Rb- 

We start with the parameterization of the Zbb interac- 
tion in Eq. ^ in terms of the shifts 8gLb,m- According 
to our fits of the electroweak data in the previous section, 
two viable regions of parameter space exist (see Table [TT] 
and Fig.[l]): 

I. 5g Rb ~ 0.016 ± 0.005, Sg Lb ~ 0.001 ± 0.001, (16) 
II. Sg Rb ~ -0.17 ± 0.005, 5g Lb ~ 0.001 ± 0.001. (17) 

For the simplest case in which the bottom quark b mixes 
with a single vector-like quark B, the shifts in the cou- 
plings are given by Sg Lb = (t 3 + l/2)s| and Sg Rb = t 3 s 2 R 
(14] . where t$ is the SU(2)l diagonal generator for the 
vector-like quark and s_L,i?. = smB^R are the sine of the 
left and right mixing angles in the b — B sector. There- 
fore, to obtain the dominant shift in 8gR b , the vector-like 



quarks must have nontrivial SU(2)l quantum numbers. 
Combined with the requirement of a charge Qb = — 1/3 
component, there are just three possible representations. 
One representation, "J ~ (3,3,2/3), allows for mixing 
with 6l and thus leads to a larger shift in Sg^b compared 
to SgRb, in contrast to the pattern indicated by the data 
( 16|17 ). Thus, there are only two viable representations: 
T 



(T,B)r 

(B,X) 



(3,2,1/6), 
' (3,2,-5/6). 



(18) 
(19) 



We now briefly review the phenomenology of each model 
in its minimal implementation. 



The first representation ( 18 ) has the same quantum 
numbers as the SU(2)l quark doublet of the SM, and 
thus contains a charge 2/3 component T and charge - 
1/3 component B. Since £3 = —1/2, mixing between b 
and B leads to a negative shift in the Zb^bn coupling, 



Sg 



Rb 



i/2- 



To resolve the A b FB discrepancy, this 



should be fixed to SgRb ~ —0.17 (17), implying a sizable 
mixing angle sr ~ 0.58. The mixing originates from the 
Yukawa couplings. The relevant mass terms in the b — B 
sector are 

-CD yiQHb R + y 2 V L Hb R + M^ L ^ R + h.c. . (20) 

The required large mixing angle 9 R is obtained if the 
Yukawa coupling is order one, Y 2 ~ 0.7M, with Yi = 
Uiv/y/2, v = 246 GeV. A similar large mixing in the 
t — T sector would induce a large Wb R t R coupling, which 
is constrained from b — > sj [2] . Thus a large breaking of 
custodial symmetry is required, inducing a sizable mass 
splitting of the T, B partners and hence a large positive 
NP contribution to the T parameter. The electroweak 
data then point to a region with light vector-like quarks 
in the 100-200 GeV range and a heavy Higgs, far above 
125 GeV. Given the observation of the 125 GeV Higgs- 
like state, as well as stringent constraints on new quarks 
in the few hundred GeV range, this model now seems 
disfavored. 



The second representation (19) with hypercharge -5/6 
contains a charge -1/3 component B, as well as an exotic 
partner X with charge -4/3. In this case t% = +1/2, so 
that the shift in the coupling 5g R b = s r/2 is positive 
and is therefore fixed to the value 5g R b ~ 0.016 (16) 



to correct A b FB . This implies a smaller required mixing 
s R ~ 0.2. Analogously to the (3,2,1/6) model, the mass 
terms in the b — B sector that lead to such mixing are 

-CD Vl QHb R + y 2 y L H^b R + M^ L ^R + h.c. . (21) 

In this simple implementation, there is no other fermion 
that can mix with the exotic quark X. One may obtain 
the required mixing angle with a relatively small Yukawa 
coupling 7/2 so that the custodial breaking is small. Thus, 
NP contributions to S and T are generically small and the 
model has an overall good agreement with electroweak 
precision data. 
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Models with vector-like fermions that couple to 
the Higgs can lead to modifications of the ef- 
fective Higgs couplings; for recent studies see 
Refs. [55], [35], @D], 03], 02, 03] Q In the simple models 
discussed here with vector-like B quarks, there is a 
reduction of the hbb coupling (by a factor c 2 R ) and to NP 
contributions to h — > gg and h — > 77 coming from loops 
of the new mirror quarks. The Higgs phenomenology for 
the (3,2,1/6) model in Eq. (20) was studied in detail 
in Ref. [IB]. In particular, for a 125 GeV Higgs the 
model predicts a large enhancement in the WW and ZZ 
channels, /Uyy ~ 2.4, in conflict with the measurements 
for the LHC and Tevatron experiments, fiyv ~ !• O n 



leading to the fermions masses is given by 



the other hand, the model in Eq. (21 ) leads only to small 



deviations from the SM, of order 10% or less 



The models described by the Lagrangians ( 20 ) and ( 21 ) 



are the minimal ones that are able to cause a sizable 
shift in the coupling Sg R b. However, the fit to the pre- 
cision data (see Fig. [I] and Eqs. ( 16|17 >), also suggests a 
small shift in SgLb , which cannot occur in the minimal 
Lagrangians. In order to bring about a small shift SgLb, 
Ref. [2] added a SU(2)l singlet vector-like fermion with 
quantum numbers (3,1,-1/3), which has a small mixing 
with the SM b^. The addition of such fermions can lead 
to additional effects in the Higgs sector beyond those 
in the minimal models due to a possibility of a large 
coupling between vector-like fermion doublet, singlet and 
Higgs. However, such a large coupling, by itself will lead 
to a sizable T parameter, due to the large splitting among 
the components of the SU (2) l doublet fermion. This can 
be easily remedied in a model with an approximate cus- 
todial symmetry. In order to exhibit the full range of 
Higgs phenomenology, we will now explore in detail an 
example of such a 'Custodial' Beautiful Mirrors model. 



Custodial Model 

We consider an extension of the SM with the following 
vector-like fermions: 



T' 
L,R 



B' 



L.R 



X' 



L.R 



(B',X') ~ (3,2,-5/6), 

(3,1,-1/3), 

(3,1,-4/3). 



(22) 



We will use primed fields to denote gauge eigenstates 
and unprimed fields to denote mass eigenstates. The B' 
and X 1 quarks have electric charges or Qb = — 1/3 and 
Qx = —4/3, respectively. 

We will allow mixing of the B vector-like fermions with 
the b quark of the SM. The most general Lagrangian 



1 See Ref. 44 for a mechanism to alter h 
charged states. 
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C D MxW^'x + M 2 B' L B' R + 
+ yi Q' L Hb' R + y 2 Q' L HB' R 
+ y^' L Hb' R + y^' L HB' R 

+ y 5 & L H^' R 
+ y^' L HX' R 

+ y 7 X' L H^' R +h. C . . 



R 



(23) 



Without loss of generality, we have rotated away a pos- 
sible mass term MB' L B' R . We note that the B' and X' 
fermions can be embedded in a doublet representation 
of a global SU(2) R symmetry, and in the limit = y§, 
2/5 = 2/7, and 2/2=2/3 = 0, the Lagrangian of the mir- 
ror fermions is symmetric under the custodial symmetry 
SU(2)l x SU(2) R which protects the model from contri- 
butions to the T parameter. To correct A b FB and Rt, we 
will need to take y 2 and y% non-zero, but still small in 
comparison to the vector-like mass terms, so that the cus- 
todial breaking is small. For simplicity, we will only con- 
sider real couplings, but see Refs. 02], 03] for recent stud- 
ies of CP- violating phases in vector-like fermion models 
which modify Higgs couplings. 

The Lagrangian p3| leads to the mass matrices in the 
B (b',B',B') and X (X',X') sectors: 



M' R 



(Y 1 Y 2 

y 3 m 1 r 4 

V Y 5 M 2 



M', 



- 52 •<"> 



where Yi = yiv/y/2. The mass matrices are are diago- 
nalzied via the orthogonal transformations: 

M B =U£M' B W B , Mx^UlM' x Wx, (25) 

Starting from the Lagrangian parameters, we can com- 
pute the mass eigenvalues in the quark sector, the ro- 
tation matrices Ub,x, Wb,x, and the couplings of the 
physical fields. However, as we will discuss in detail mo- 
mentarily, the masses of the new quarks are quite con- 
strained by LHC searches, and as such must be heavy. 
It is therefore useful to work in the regime where the 
vector-like mass terms Mi 2 3 are much larger than the 
corresponding Yukawa terms Y i: and integrate out the 
heavy vector-like fields. We can then match to a set of 
effective operators that can be used to describe the cor- 
rections to the SM predictions for e.g. the Zbb couplings 
and Higgs couplings. We now describe these steps and 
the approximate results in the effective theory. 

Effective theory We integrate out the new vector- 
like fermion fields at tree level, in the regime of large 
Mi 2,3- The procedure is as follows: we use the classical 
equations of motion to solve for "J, B, X in terms of 
the light fields Q,bu, and substitute these solutions back 
into the Lagrangian. We then make use of the classical 
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equations of motion for Q, b R to match to the following 
operator basis: 





= i(H^D^H)(b R ^b R )+ h.c, 


(26) 


S H Q 


= i(HtD IA H)(&fQ) + h.c, 


(27) 




= t(H^ a D^H)(Q^a a Q) + h.c., 


(28) 


Ohy 


= (H^H){QHb R ) + h.c, 


(29) 


O y 


= {QHb R ) + h.c.. 


(30) 



The effective Lagrangian is given by 



(31) 



The first three operators above ( 26|27|28 l cause modifi- 
cations to the Zbb vertex, given by 



$9Lb = 



(a 



HQ 



HQ 



)v 2 



5gm = - 



a H bV 



(32) 



The last two operators equations above ( 29|30 1 affect the 
bottom quark mass and the hbb coupling: 



m b 



Uhbb 



ayv 

71" 

av 



2%/2 
Sclhyv 2 
2^/2 



(33) 
(34) 



After integrating out the heavy fermions and match- 
ing to the operator basis above, we obtain the following 
coefficients for the operators: 



am 



ay 



2Mf 



-yi, 



l HQ 



O-HY 



4Mf ' 



a HQ 



4M 2 2 ' 



Vivi , Vwl 



2M\ 



IMl 



M 1 M 2 



(35) 



We thus obtain the approximate expressions for the bot- 
tom quark mass, the hbb coupling, and the shifts to the 
Zbb vertex: 



Using Eq. (38), we can rewrite the approximate ex- 



pression for the bottom mass (36 land hbb coupling (37) 
in terms of the shifts £<7l6, Sg Rb . Since 5gLb,m <C 1, we 
can express r b , the ratio of the h —t bb partial width to 



its value in the SM, defined in Eq. ( 10 ) as 



n 



Uhbb 

m b /v, 
1 - 4(5 9m 



(40) 



SgLb) + 8^Sg Rb Sg 



where we have taken Y 2 , Y 3 to be positive. Thus, we see 
that only for Y5 3> m b can r b deviate significantly from 
unity. 

Low energy Higgs theorems Here we will apply 
the low energy Higgs theorem [45] , [46] to derive approx- 
imate expressions for the corrections to h 77 and 
h — > gg. In all generality, integrating out heavy fermions 
/ charged under electromagnetism (color) yields the fol- 
lowing effictive Lagrangian: 



C D 



167TV 



O.s 



E >t 

f=B,X 



b™Z f 



E bC f z f 

f=B,X 



1&1TV 

where the quantities Zf are given by 
d 



hGl v G^\ (41) 



z f 



dlogv 



E lo s m /,*( u ) 



(42) 



and the sum is over all heavy fermions. The beta 
function coefficients are in general b^ M = (4/3)-/V c Qj, 
b c j = (4/3)C(r), where N c , Qj and C(r) are the number 
of colors, the electric charge and the Casimir of the par- 
ticular representation of the fermion /. For the system 
at hand we have 



m b 



Yi 1 



Uhbb 



Yi 1 



V 2 

2M| " 
3Yj 
2Mf 



SgLb 



Y' 2 

2M| ' 



Y 2 
2M 2 
ZY 2 
2M 2 



Sg Rb 



Y 2 Y 3 Y 5 
MiM 2 ' 
3Y 2 Y 3 Y 5 
M\M 2 



Y 2 

2M 2 ' 



(36) 
(37) 
(38) 



Corrections to Zbb and hbb couplings As dis- 
cussed above, the mirror quark representation (3,2,-5/6) 
will lead to a strictly positive shift in the coupling 5g Rb , 



and this is seen in Eq. (38 1. Thus, we will fix the shifts 
SgLb and 5g Rb to their best-fit value in Eq. ( 16 1 in order 



to bring the A b RB and R b predictions into accord with the 
experimental values. We then obtain from Eq. (38) the 



following relations between the Lagrangian parameters: 
F 2 ~±0.04M 2 , ¥"3 ~ ±0.17 Mi. (39) 



,EM 
°B 



4/9, b x M = 64/9, b B =b c x = 2/3. (43) 
Given the effective Lagrangian, we can derive the ratios 



1 B Zb 



b x z x 



ASM 

b c B z B + b c x z x 
Asm 



(44) 



With 



1.3 and A^ M 



6.6. 



To evaluate the quantities Zf, it is convenient to rewrite 
the sum in terms of the determinant of the mass matrix, 
so that it can be evaluated in any basis. For the bot- 
tom sector, however, we cannot straightforwardly do this 
since the mass matrix contains a light eigenvalue. Let 
us then derive a simple modification which allows us to 



apply the low energy theorem. We can rewrite Eq.(42) 
as 



ZD 



d 



3 



d log v ^— ' 

b i=2 



9 

log det M B M B 



2v d 



m b (v), (45) 



dlogv ' rrib(v) dv 

where the sum on the first line is only over the heavy 



eigenvalues of the mass matrix M.' defined in (24) and 



is the usual piece associated to the use of the low energy 
theorem. Since the determinant includes a light bottom 
quark contribution, we must subtract this contribution, 
which is accounted for by the second term above. Using 



the approximate expression for mb in Eq. ( 36 1 , we can 
write the leading terms of order 1/M 2 : 



zb 



2Y£ 2Y 3 2 



M| 



M 2 M X M 2 



4Y 4 Y 5 



Ml M 2 ™ 

Note that the last expression is valid in the regime 
1*4,5 3* ii,2- In a n y case, since the first two terms above 
are directly proportional to the small shifts in the Zbb 



couplings ( 38 ) , they will not cause large deviations in 
the r gn . 

For the X sector, we can straightforwardly convert the 
sum to a determinant. The result is 

Q 

log det MxMx 



(47) 



Zx d log v 
~ M x M-s ' 

Note that, due to its large electric charge of —4/3 (lead- 
ing to b^ M — 64/9), the X fermion contribution will be 
the dominant one in h 77, and correspondingly Ys,7 
will mainly affect r 7 . On the other hand, the B and X 
contributions to r g will be comparable. 

Oblique corrections The new mirror fermions in 
our model can give sizable contributions to the oblique 
parameters S and T, especially if they have large Yukawa 
couplings to the Higgs field. Given that we are invoking 
NP shifts to the couplings SgLb,Kb to correct A b FB and R b , 
the allowed regions for S and T will be shifted somewhat 
compared to the usual regions derived for the SM Zbb 
couplings. 

We have thus performed a fit to the electroweak data 
allowing for the presence of coupling shifts 5 g Lb, Kb and 
oblique parameters S and T. For the Zbb coupling shifts, 
the fit selects regions essentially identical to the one dis- 
played in Fig. Q, while for the oblique parameters the 
preferred regions are 



S = -0.02 ±0.09, T = 0.03 ±0.08, 



(48) 



with a correlation coefficient of 0.90. Since the motiva- 
tion for introducing the mirror fermions is to improve the 



description of the precision electroweak data, we will re- 



strict the parameter space to the region in Eq. ( 48 1 when 
we present our numerical results below. 

For the model under consideration, the oblique cor- 
rections are most easily computed directly in the mass 
basis. Our procedure will thus be to diagonalize the sys- 



tem ( 23 ) , obtain the physical couplings of mirror fermions 



to gauge bosons, and compute the vacuum polarization 
functions of the gauge bosons. Since the physical cou- 
plings of the top and bottom quarks to the electroweak 
gauge bosons are shifted in this model from the SM, 
we will include in our computations the contributions of 
those particles and at the end subtract off the SM values. 

Here we present the general expressions for the vacuum 
polarization functions. Given two fermions ipiy^Pj inter- 
acting with a vector boson V A through the Lagrangian 



(49) 



the contribution to the vacuum polarization function 
H l ^ g (p 2 ) arising from the one-loop exchange of these two 
fermions is given by: 



AB ~ 16tt2 



(50) 



where wc have defined 



AB 



r'JrJ< 1 d*J dJ* 
^A B + n A n B ' 



F(p,m l ,m J ) = 4:B 2 2{p,mi,mj) + 2p 2 (B 2 i(p,m l ,m :j ) 



B 1 (p^ mj ))^mj- m ] + "-, 



G{p,m l ,m :j ) = -2B (p, m u m 2 ), 



(51) 



with Bi(p,mi,m,j) are the Passarino-Veltmann func- 
tions |47j . With the vacuum polarization functions in 
hand, we can straightforwardly compute S, T according 
to the definitions in e.g. [15]. 

Collider bounds and prospects The collider phe- 
nomenology of this model is very similar to the 'mini- 
mal' model with the (3, 2, —5/6) representation discussed 
above, and has been studied in detail in Ref. [T7]. For 
mB,x 700 GeV, QCD pair production has the largest 
rate, while for larger masses electroweak single produc- 
tion of mirror quarks dominates. 

As recently summarized in |49j , numerous searches for 
heavy vector-like quarks have been performed both at the 
Tevatron and LHC, though no evidence for the existence 
of exotic quarks has been observed. Here we focus on the 
constraints and signatures relevant for our model. 

The exotic bottom-type quarks B arising in our model 
will be pair-produced at hadron colliders via their QCD 
interactions. Since they only mix with the third gen- 
eration down-quark of the SM, their decay modes are 
B -> tW, bZ, bh. 
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The first decay mode B — > tW is the most constrained 
by the experimental searches for vector-like 6-quarks. 
Searches have been performed with final states contain- 
ing only one [501 EE] (two [H [53] or at least two [MUSS]) 
W bosons decaying leptonically. The strongest constraint 
comes from the ATLAS analysis [55] based on 4.7 fb _1 
7 TeV LHC data. This analysis requires final states con- 
taining at least two isolated leptons of the same charge 
and at least two jets, including at least one b-tagged jet: 
vector-like 6-quarks decaying 100% to tW with mass be- 
low 670 GeV are now excluded at 95% C.L. . 

Furthermore, several Tevatron and LHC searches give 
bounds on the mass of vector-like 6-quarks decaying 100% 
to bZ. An early CDF analysis based on 1.06 fb _1 56 of 
data bounds the mass to be greater than 268 GeV. This 
bound has recently been improved by the ATLAS [57J 
and CMS collaborations [58], using respectively 2 and 
4.9 ftr 1 data at the 7 TeV LHC. Vector-like 6-quarks 
decaying 100% to bZ with mass below 550 GeV are now 
excluded at 95% C.L. 

Finally £?-quarks decaying to bh, giving 66 final states 
are the most difficult to search for. The ATLAS and 
CMS collaborations are performing extensive searches for 
stops and gluinos. Several analysis, based on 7 and 8 TeV 
LHC data, are requiring final states containing multiple 
jets and at least 1,2,3 6-jets [59"H6"3"] . This searches can be 
recast to set bounds on the mass of vector-like i?-quarks 
decaying mainly to bh. However, to optimize stop and 
gluino searches, the CMS and ATLAS analysis impose 
a strong requirement on the amount of missing energy 
(at least 100-150 GeV), which supresses the sensitivity 
of these searches to £?-quarks decaying to bh. The con- 
straints on the bh decay mode are thus weak in compar- 
ison to the tW, bZ modes, thus warranting a dedicated 
6b search. 

The additional exotic X quarks arising in our model 
will also be produced copiously in pairs and will decay to 
bW~ and BW~ . Tevatron and LHC t' searches thus set 
stringent bounds on the mass of these vector-like fermions 
decaying to bW . In particular, searches have been per- 
formed with one [641466] and two [67] leptonic W bosons. 
The CMS analysis based on 4.6-4.7 fb" 1 7 TeV LHC 
data [5^ sets the most stringent bound: M x > 560 GeV. 

The longer decay chains X BW -> (tW, bZ, bh)W 
are instead more weakly constrained by present Tevatron 
and LHC analysis. However, it is quite difficult in our 
model to arrange for the decay X — > BW. As we will 
see below, the Higgs signal strength data prefer Yukawa 
couplings to the X sector that are larger than those to 
the B sector, typically causing the lightest mirror quark 
to be X rather than B. One could arrange a lighter X by 
making A/3 <C M2, but at the expense of large custodial 
symmetry breaking. Thus in practice, the exotic X quark 
in our model decays via X —> bW with a 100% branching 
ratio and the constraint from t' searches, Mx > 560 GeV 
is quite robust. 



Given that X quarks must be heavier than 560 GeV 
and assuming Mx < Mb, as indirectly suggested by 
Higgs signal strength data and custodial symmetry, in 
practice the bounds on Mb from vector like 6-quark 
searches are automatically circumvented in our model. 
In particular the bound Mb > 670 GeV obtained by 
CMS for 6-quarks decaying 100% in tW is satisfied since 
the B quarks typically have 0(1) branching ratios to bZ 
and bh, weakening the limit. 

In terms of future prospects, we estimate that with the 
full 8 TeV data set, the bound on X coming from pair 
production can be extended to roughly mx > 700 GeV. 
To probe higher masses, it will be important to focus on 
single production of mirror quarks. We expect that the 
= 14 TeV LHC will be able to probe mirror quarks 
masses into the ~ 2 TeV range; see [T7J for a detailed 
study. 

Finally, we mention that besides the direct production 
of vector-like quarks, another possible novel signature is 
enhanced double Higgs production, which occurs due to 
the loop-level contribution of the heavy quarks to gg —¥ 
hh. See, for example, Refs. [55] for dedicated studies of 
this signature. 

Vacuum Stability The Yukawa couplings of the 
mirror fermions to the Higgs boson yield negative con- 
tributions to the /3 function of the Higgs quartic cou- 
pling. For O(l) Yukawa couplings, this can cause the 
Higgs quartic to run negative at a low scale, potentially 
leading to a vacuum instability problem. This point has 
been emphasized recently [39] , [40] , [41] in the context of 
the enhanced the di-photon signal. In the model under 
consideration this is also an issue as can be seen by ex- 
amining the (3 function for the Higgs quartic A: 

24A 2 + 

-6(y t 4 



16tt 2 /3 a 



12A(y t 2 



vl 



vl 



y 2 7 ) 



(52) 



This expression is valid in the limit of 2/1,2,3 <C yt and of 
small gauge couplings. On the one hand, 0(1) values for 
the new coupligs 1/4,5,6,7 can cause desirable modifications 
to Higgs physics, but clearly the larger these couplings 
are the faster the Higgs quartic will run negative. As an 
example, if we take the values 2/4 = 2/5 = and 2/6 — 
2/7 = 1 and a threshold at the vector-like mass scale of 
M = 800 GeV, we find that the Higgs quartic vanishes at 
the scale 2 TeV, while it can be metastable up to scales 
of order 3 TeV, where we have used the results of [69 to 
estimate the metastability bound. 

Thus, the model with these parameters requires a UV 
completion at low scale. One possible completion would 
be supersymmetric version of the model, since the scalar 
superpartners of the mirror quarks can provide an equal 
and opposite contribution to the beta functions. This 
direction, however, lies outside the scope of this paper. 

Numerical Results We now explore the parameter 



space of the model ( 23 ) , finding regions which resolve 
A b FB and Rb discrepancies, give small contributions to 
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the oblique parameters S and T, are consistent with di- 
rect searches for vector-like quarks at Tevatron and LHC, 
and, finally, provide a good description of the Higgs sig- 
nal strength data. 

Since there are many new parameters in the model, we 
will now make several physically motivated assumptions 
in order to reduce the parameter space: 1) We fix the 



couplings Y 2 and Y 3 according to Eq. (39) in order to 



cause t he s hift in the Zbb couplings to their central values 



in Eq. ( 16 1, which bring the A FB and Rb predictions into 



agreement with their measured values. 2) We fix Y\ by 
the requirement of obtaining the correct b quark mass via 



Eq. (36) 3) We fix Af 2 = M 3 as motivated by custodial 
symmetry. This still leaves 5 parameters that describe 
the model, namely, Y4,Y 5 ,Y e ,Y T , Mi, M 2 . For simplicity, 
we will further fix I4 = Y 5 = Yb and Y e — Y 7 = Yx, 
describing common Yukawa couplings in the mirror B- 
quark and V-quark sectors, respectively. Finally, we will 
also assume a common vector-like quark mass scale Mi = 
AI2 = Ms = M. With these simplifications there are 3 
parameters: Yb,Yx, and M. 

The most robust constraint on the model comes from 
the t' searches at the LHC discussed above, which re- 
strict Mx > 560 GeV. With the simplifying assumptions 
above, the lightest X mass is given by 



M X =M - Y : 



x- 



(53) 



so that, e.g. for a given vector-like mass scale M, the X 
sector Yukawa coupling must be less than some maximum 
value. 

Wc first fix the vector-like mass scale to M — 800 GeV 
and show in Fig. [3] the preferred regions of parameter 
space in the Yb — Yx plane. The fit to the Higgs signal 
strength selects a region -100 GeV < Y B < 20 GeV 
and \Yx\ > 100 GeV, shown in dark blue in Fig. [3j In 
this region, the signal strength in the diphoton channel is 
enhanced, 1 < /i 77 < 1.6, as illustrated by the constant 
/i 77 contours in orange in Fig. [3] This enhancement in 
result of 1) the loop contribution of the charged 
—4/3 particle X to h — > 77 causing r 7 > 1, and 2) a 
suppressed coupling of the physical 6-quark to the Higgs 
due to mixing, which causes < 1 and enhances the 
branching ratio of h — > 77. The brown shaded region 
in Fig. [3] corresponds to Mx > 560 GeV and is thus 
excluded by t' searches. Finally, the parameters in the 
gray shaded regions lead to oblique contributions which 
are in tension with the preferred values of our fit, given 



in Eq. (|48). 

Next, we fix the common £?-sector Yukawa coupling 
to the value Yb = —65 GeV and show in Fig. [4] the 
preferred regions of parameter space in the M—Yx plane. 
The regions are as described above for Fig. [3] We see 
that, for this value of Y B , one obtains an enhancement of 
/i 77 ~ 1.5, though only a small region of parameter space 
is allowed. Decreasing Y B , larger regions of parameter 
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FIG. 3. Preferred regions of parameter space in the Yb — Yx 
plane, with the common vector-like mass fixed to M = 800 
GeV. The dark (light) blue area represents the Higgs signal 
strength la (2a) preferred region. The gray region predicts 
oblique parameters 5*, T outside the la preferred region. The 
brown shaded region is excluded by t' searches. We also 
show in orange the contours of constant signal strength in 
the diphoton channel /x 77 . 



space open up, although the enhancement in ^ 77 is less. 
This is consistent with the regions displayed in Fig. ([3|. 

For these results, we have taken V 2 and Y3 to be posi- 
tive, so that the Higgs data selects out regions of negative 
Yb which suppress r\, and thus enhance /i 77 . The depen- 
dence of rb on the Yukawa couplings can be seen in in 
Eqs. ( 36|37|40 ). We note that there are several physi- 
cally equivalent regions related by different sign choices 
of the Yukawa couplings. 

We concluded that the model can provide a good de- 
scription to the Higgs signal strength data, and in partic- 
ular an enhancement in the diphoton channel /i 77 ~ 1.6, 
while simultaneously resolving the two discrepant preci- 
sion observables A b FB and Rb- 



CONCLUSIONS 

In this paper we have reexamined the precision elec- 
troweak data following the discovery of a Higgs-like state 
at the LHC. We have paid special attention to A b FB and 
Rb, which currently exhibit tension with the SM at the 
level of 2.4<t and 2.3cr, respectively . While the A b FB 
discrepancy is well-known, the tension in Rb is a recent 
development due to a state-of-the-art calculation includ- 
ing two-loop electroweak corrections [H]. 
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FIG. 4. Preferred regions of parameter space in the M — Yx 
plane, with the common B-sector Yukawa fixed to Ya = —65 
GeV. The dark (light) blue area represents the Higgs signal 
strength la (2a) preferred region. The gray region predict 
oblique parameters S, T outside the la preferred region. The 
brown shaded region is excluded by t' searches. We also show 
in orange contours of constant signal strength in the diphoton 
channel. n 11 



One logical possibility is that these discrepancies are 
the result of an unknown systematic error, although no 
candidate for such an effect is known. Nevertheless, we 
have shown that upon exclusion of these observables from 
the global electroweak fit, the remaining data are well 
described by the SM, and the fit prefers a Higgs mass of 
125.7 GeV due to the precise LHC measurement. While 
no strong argument for new physics can be made in this 
case, a marginal improvement in the fit is possible if NP 
causes oblique corrections. 

On the other hand, if one believes the measurements 
are correct, then the global fit to the SM is of low quality. 
New physics in the form of non-universal shifts to the Zbb 
vertex can dramatically improve the global fit and bring 
the predictions for A b FB and into agreement with the 
measurements . 

It is reasonable to expect that such new physics will 
lead also to modifications to the properties of the Higgs 
boson, and this interplay has been the primary focus of 
our work. In particular, we have investigated the Beau- 
tiful Mirrors scenario, which contains new vector-like 
quarks that mix with the bottom quark of the SM. This 
scenario generically predicts modifications to the Higgs 
boson couplings to b quarks, gluons, and photons. We 
have derived in a model-independent fashion the shifts 
preferred by the Higgs signal strength data from the LHC 
and Tevatron experiments. 



In the simplest models in this framework, the new mir- 
ror quarks in this scenario can have large couplings to 
the Higgs boson, but only at the expense of large custo- 
dial symmetry breaking which could reintroduce a ten- 
sion with the electroweak data. With this motivation, we 
have proposed a 'custodial' version of the model which 
can protect against large contributions to the T param- 
eter. We have confronted this model with precision elec- 
troweak data, Higgs signal strength data, and collider 
searches, deriving regions of parameter space which are 
in agreement with all experimental results. In particular, 
the model predicts an enhancement of the diphoton sig- 
nal strength between 1 < /x 77 < 1.6. If the current trend 
in the diphoton channel continues, then the scenario is 
testable in the immediate future, as it predicts light mir- 
ror quarks around the TeV scale which will be probed by 
the LHC experiments. 
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